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We examine measurements of time-dependent CP asymmetries that could be made in new and 
future flavour facilities. In charm decays, where they can provide a unique insight into the flavor 
changing structure of the Standard Model, we examine a number of decays to CP eigenstates and 
describe a framework that can be used to interpret the measurements. Such measurements can 
provide a precise determination of the charm mixing phase, as well as constraints on the Standard 
Model description of CP violation an possible new physics contributions. We make a preliminary 
assessment, based on statistical considerations, of the relative capabilities of LHCb with data from 
pp collisions, with Belle II and SuperB using data from Bd, B s and charm thresholds. We discuss 
the measurements required to perform direct and indirect tests of the charm unitarity triangle 
and its relationship with the usual Bd triangle. We find that, while theoretical and experimental 
systematic uncertainties may limit their interpretation, useful information on the unknown charm 
mixing phase, and on the possible existence of new physics can be obtained. We point out that, 
for Bd decays, current experimental bounds on Ars d will translate into a significant systematic 
uncertainty on future measurements of sin 2/3 from b —¥ ccs decays. The possibilities for simplified 
B s decay asymmetry measurements at SuperB and Belle II are also reviewed. 
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I. INTRODUCTION 

The Standard Model (SM) description of quark mixing 
and CP violation is described by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [TJ[2]- This matrix can be writ- 
ten as 

(v ud v us v ub \ 

Vckm = V cd V cs V cb , (1) 

V Vta v ts v tb J 

where the Vij are coupling strengths for up-type to down 
type quark transitions. Unitarity of the CKM matrix 
gives rise to six triangles in a complex plane, one of which, 
the bd triangle, has been extensively studied by the B 
Factories, and has earned the name of 'The unitarity tri- 
angle'. The unitarity triangle is 1 

V* d Vub + V c * d V cb + V t * d V tb = 0. (2) 

The angles of the unitarity triangle are a = (91.4±6.1)°, 
/3 = (21.1 ± 0.9)°, and 7 = (74 ± 11)° g]. Given 
that, for three generations, these angles must add up 
to 180°, the measurements of a and /3 provide a strong 
constraint on the value of 7 in the SM. Precision tests 
of the CKM mechanism have been made only for tran- 
sitions of down type-quarks from the second and third 
generations, and time- dependent CP asymmetries have 
only been measured in the third generation (B decays). 



In order to study the corresponding phenomena with an 
up-type quark one has to study the charm decays as top 
quarks hadronize before being able to form a quasi-stable 
meson. In this paper we concentrate on examining ap- 
proaches that may be usable, at existing and future ex- 
perimental facilities, in order to study time-dependent 
CP asymmetries in the charm sector, and how results 
might be interpreted in the context of the CKM matrix. 
In particular, we focus on ways in which the SM expec- 
tations of the "charm triangle" , defined in Eq. ^ might 
be examined — something that has yet to be done. In 
addition to this, we make a few observations on measure- 
ments of B d and B s decays in Sections |VIII| and |IX| 

In addition to the bd triangle, unitarity of the CKM 
matrix also gives rise to the charm (cm) triangle 

v: d v cd + v: s v cs + v* b v cb = 0, (3) 

which depends on the weak phase 7 by virtue of the 
presence of the factor V ub . The angles of the charm 
triangle can be written as a c , /3 C , and 7 C . Some time 
ago Bigi and Sanda [5] pointed out that 7 C ~ 7, and 
a c = 180 o ~7 c + O(A 4 Y= 180° -7 + C(A 4 ) 2 . Hence the 
charm and unitarity triangles are related in a simple way. 
Using the Buras et al. variant of the Wolfenstein param- 
eterization [7] of the CKM matrix up to 0(A 5 ), the 
weak phase f3 c can be estimated to be ~ 0.035°. Hence 
the sum of a c and 7 C should be essentially 180°. Existing 



2 Our nomenclature differs from that used in Ref. Our angles 
We depart from the usual convention by defining the complex can be related to theirs as follows: fi c = (j>g u , -f c = <I>2 U > an d 

conjugates of the triangle sides. a c = (p'p' . 
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constraints on the Wolfenstein parameters can be used to 
give a clean prediction of the cu triangle parameters. In 
order to verify if the CKM matrix is the correct descrip- 
tion of quark mixing the angles a c , f3 c , and 7 C need to be 
measured, as well as the sides of this triangle. The e + e~ 
collider experiment SuperB is the only facility where one 
can, in principle, perform all of the necessary measure- 
ments to perform a complete cross-check of the two trian- 
gles. This requires large samples of B, D and D s mesons, 
which Superi? will accumulate through runs at charm 
threshold, and at the T(45). In order to interpret time- 
dependent measurements in terms of angles of the cu tri- 
angle, a precise measurement of the charm mixing phase 
is required. We propose that one studies D — > K + K~ to 
measure this mixing phase, and the difference between 
measurements of D — > K + K~ and D — >• tt + tt~ will then 
give —2j3 Cie ff. A single measurement of, or constraint on, 
P c , predicted to be (0.0350 ± 0.0001)°, would clearly be 
of interest, but this will require a careful study of effects 
of other amplitudes and possible long range effects. Any 
observed deviation from this expectation would then be 
an indication of new physics (NP). Indeed it is worth not- 
ing that the measurement of sin 2/3 from B meson decays 



to final states including Charmonium and a neutral kaon 
are inconsistent with the SM at the level of 3.2a [8 . This 
result is a strong motivation to perform the correspond- 
ing studies of the SM in the charm sector, which is the 
subject of this paper. 

The possibility of large CP violation effects in charm 
decays has been discussed elsewhere [9Tll 2], however un- 
til now these have focused on time-integrated measure- 
ments, and ignored possible time-dependent effects. It 
is clear that we need precision experimental tests of the 
unitarity triangle, in particular the angle 7, and the sides 
of both the charm and unitarity triangles. In addition to 
this we also need to start measuring the charm triangle 
angles precisely in order to validate the CKM descrip- 
tion of CP violation for up-type quarks. Theoretical un- 
certainties will ultimately limit the constraints that can 
be placed on the SM, and we discuss some of the issues 
here. The remainder of this paper outlines the details re- 
quired to perform time-dependent CP measurements in 
the charm sector, using CP eigenstate decays, and con- 
straints on the sides of the charm triangle, before making 
a few observations on B d and B s decays. 



II. THE CKM MATRIX 



The CKM matrix given in Eq. ([I]) can be parameterized in a number of different ways. The Wolfenstein parame- 
terization [B] is an expansion in terms of A = sin# c , A, p, and i], where 9 C is the Cabibbo angle. A variant on this 
parameterization has been proposed Buras et al. 7], and has the advantage of preserving unitarity to all orders in A 
in the "fed" triangle, though possibly not in the "cu" triangle. The Buras et al. variant of the CKM matrix up to and 
including terms C(A 5 ) is given by 

/ l-A 2 /2-A 4 /8 A A\ 3 (p-n 1 )\ 

Vckm= [ -X + A 2 X 5 [l-2(p + ir 1 )}/2 l-A 2 /2-A 4 (l + 4A 2 )/8 AX 2 \+0(X 6 ). (4) 

\AX 3 [l-(l-X 2 /2){p + i V )} -AX 2 +AX i [l-2{p + ir 1 )}/2 l-A 2 X i /2 ) 

The choice of convention used to interpret data in terms of physical observables is irrelevant as long as sufficient terms 
in the expansion are used. Expansions to C(A 3 ), have been sufficient for the B factories era, however one should 
consider additional terms as we move into the era of LHCb and the Super Flavor Factories (SFF's). 

The apex of the unitarity triangle, obtained when Eq. ^ is normalized by V cd V* b , is given by the coordinates (p, rj), 
where 



p = p[l-A 2 /2 + 0(A 4 )], (5) 

fj = V [1 - A 2 /2 + 0(A 4 )] . (6) 
The CKM matrix may be written in terms of p and fj as 

l-A 2 /2-A 4 /8 A AX 3 {p-ifj){l + X 2 /2) 

Vckm = ( -A + A 2 X 5 [l-2{p + if})}/2 1- A 2 /2- A 4 (l + 4A 2 )/8 AX 2 \+0(X 6 ). (7) 

AX 3 [1 - p - ifj] -AX 2 +AX 4 [l-2{p + ifj)}/2 l-A 2 X*/2 



Current constraints on the CKM parameters A, A, p, rj, 
p, and fj from global fits [3l [13] are given in Table |TJ 

The angles of the unitarity triangle given in Eq. (§ are 



a, P, and 7, where 

a = arg [-VM/VM] = (91.4 ± 6.1)°, (8) 

P = arg {-V cd V; b /V td V t i} = (21.1 ± 0.9)°, (9) 

7 = arg [-V ud V* b /V cd V* b ] = (74 ± 11)°. (10) 
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TABLE I: Constraints on the CKM parameters A, A, p, r\, 
and fj obtained by the UTFit and CKM fitter groups. 



(a) 



Parameter 



UTFit 



CKM Fitter 



Mean Used 



A 
A 
P 
n 
p 



0.22545 ± 0.00065 0.22543 ± 0.00077 0.22544 ± 0.00705 



0.8095 ± 0.0095 
0.135 ±0.021 
0.367 ±0.013 
0.132 ±0.020 
0.358 ±0.012 



0.812 



+ 0.013 
0.027 



0.144 ±0.025 
■1A2+ 016 



0.811 ±0.015 



0.138 ±0.022 
0.350 ±0.014 



v* d v ub 



v t * d v tb 




v; d Vc b 



The most precisely measured angles are a and f5 using 
B meson decays into pp |141 115j and charmonium final 
states [HJ [17] , respectively. Given unitarity, in the SM 
with just three generations, only two of these angles are 
independent, hence 7 is, in principle, a redundant cross- 
check of the CKM matrix. 

These angles can also be computed from values and 
uncertainties for A, A, p and fj. Taking simple averages 
of CKM Fitter and UTFit values in Table [EJ the angles, 
computed to order A 6 , are 



(b) 




V* V 

v us v cs 



v: b v cb 



FIG. 1: (top) The bd unitarity triangle in Eq. Q, and (bot- 
tom) the cu unitarity triangle of Eq. Q. 



a = (89.4 ±4.3)°, 
P = (22.1 ±0.6)°, 
7 = (68.4 ±3.7)°. 



(11) 
(12) 
(13) 



Comparing Eqns ^ and ( fl"0| ) with Eq. Q, one can 



sec that V* 



td 



\Vtd\e 



-03 



and V ub 



V ub e- 



These 



relations are exact for low orders of A, and the equality 
breaks down as V cd is complex at order A 5 . 

The angles of the charm unitarity triangle given in 
Eq. (|3| are 



7c 



arg| 
arg| 
arg| 



-v: b v cb /v: s v cs \ 
-v: d v cd /v: s v cs \ 
-v: b v cb /v: d v cd ] 



(14) 
(15) 
(16) 



where as already noted j c ~ 7 and a c = 180° — 7 ± 
C(A 4 ). Again, using the averages of CKM Fitter and 
UTFit values for A, A, p and fj and their errors, we predict 
that, to order A 6 



a c = (111.5 ±4.2)°, 

P c = (0.0350 ±0.0001)°, 

7c = (68.4 ±0.1)°. 



(17) 
(18) 
(19) 



These predictions for the angles of the charm triangle 
could, and should, be tested experimentally, either di- 
rectly (through time-dependent CP asymmetries) or in- 
directly (through measurements of the sides of the trian- 
gle). On comparing Eq. (151 with Eq. Q, one can see 
that V cd = \V cd \e l (P c ~~ 7T \ Both the bd and cu triangles 
are shown in Fig. [l] 



III. TIME-DEPENDENT EVOLUTION 

Neutral meson mixing is a phenomenon that only oc- 
curs for K, D, and B d s mesons (Charge conjugation is 
implied throughout). Here, in describing the formalism 
common to these systems, we refer to the mesons as P. 
The effective Hamiltonian describing neutral meson mix- 
ing is given by 



= M- 



eff 



Hence neutral meson mixing can be described by 
. d 



f M11 


M i2 A 




Tu r 12 \ 


1 M31 


M 22 ) 




v ^21 ^22 J 



\p°) 

^0, 



i dt\\p*)) = r-2 r Mip 



p°) 
0, 



(20) 
(21) 

(22) 



where \P°) and \P ) are strong eigenstates of neutral 
B, D, or K mesons. The matrix elements in Eq. (22) 
must satisfy Mi% — M 22 and Tn = T 22 in order to be 
consistent with CPT symmetry. A further constraint can 
be obtained in the limit of CP or T invariance, where 
ri 2 /Mi 2 = r 2 i/Af 2 i must be a real quantity. 

One can write the mass eigenstates as an admixture of 
the strong eigenstates in the following way 



\P 1 , 2 )=p\P*)±q\P"), 



(23) 



where q 2 ± p 2 = 1 to normalize the wave function, and 



q 
p 



1 M- 



12 



iT* 12 /2 



M 12 - iT 12 /2 



(24) 
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The magnitude of q/p is very nearly one in the SM. If one 
considers the mass eigenstates under the CP operator, it 
follows that |Px) is CP even, and |P 2 ) is CP odd. The 
mass and width differences AM and Ar between the 
mass eigenstates are given by 



where neutral mesons oscillate from particle to anti- 
particle state with the characteristic mixing frequency 
AM. Detailed discussions of this formalism can be found 
in a number of text books. 



AM 

Ar 



M 2 -Mi, 

ri-r 2 , 



(25) 
(26) 



A. Uncorrelated meson production 

It can be shown that the general form of the time-evolution of a neutral meson decaying into some final state / is 
given by 



T(P° -»•/)« e - rit 



1 + e 



Am 



T(P U -> /) cx e 



-Tit 



2 



1-e 



aj ■ / ^ e Art/2 I cog AM ^ _ 2Jmgy s . n 



1 + |A/I 2 



2 + l + |A,| a ^ j >v 1 + |A,| 2 +1 + IA,! 2 smIXMt J 



,(27) 
,(28) 



where 



zi 4 



(29) 



and A (A) is the amplitude for the P (P) decay to a 
final state /. Note that X f is not related to the CKM ex- 
pansion parameter A discussed above, but is a complex 
parameter related to mixing and decay transitions. The 
time t = is defined by the production of a definite me- 
son state (flavor, CP or mixed flavor), that subsequently 
evolves as a P — P admixture until it too decays. The 
identification of the flavor of a meson state at some fixed 
point in time is critical for a time-dependent measure- 

" if \q/p\ ± i 



then 
A 1 2 there 



ment and is discussed in Section IV 
there is CP violation in mixing, and if \A\ 2 ^ 
is direct CP violation, hence a measurement of the real 
and imaginary parts of Xf (or equivalently the magnitude 
and phase) is able to probe the combination of these two 
effects i.e. interference between mixing and decay. It 
should be noted that for all time-dependent CP asym- 
metry measurements of decays made by experiments 
until now the assumption that Ar = has been used. 
We discuss the ramifications of this assumption in Sec- 
tion Ell 

A time-dependent decay rate asymmetry can be com- 



puted from Eqns ( 27 ) and ( 28 ) as follows 



A(t) 



nt)-T(t) 
r(t) + r(ty 



2e ATt/2 



(|A, 



1) cos AMt 



(30) 

2ImX f sin AMt 



(1 + lA.PXl + e Art ) + 2PeA / (l - e Art ) ' 



where T(t) and T(t) are the time-dependent rates, re- 
spectively, for P° — >• / and P° —> f transitions. The 



asymmetry depends on the real and imaginary parts of 
X f as well as |A f | 2 , hence it is possible to extract X f from 
data in terms of only two parameters, the real and imag- 
inary parts of X f , as |A f | 2 is completely correlated with 
those parameters. This formalism is normally written in 
terms of hyperbolic functions, and one can derive those 
results by combining the exponential factors in the equa- 
tions above. In the limit that Ar = 0, Eq. (30) reduces 
to the familiar result 



A(t) = -C cos AMt + S sin AMt, 



where 



S 



2ImX f 



and C 



1+ A f 



(31) 



(32) 



This approximation has been used in the B factory mea- 
surements of the angles in the bd unitarity triangle. For 
future measurements, we note, the validity of the as- 
sumption that Ar = will need further checking. 



B. Correlated production of neutral mesons 

Neutral K, D, or B mesons are produced in correlated 
pairs in e + e~ collections with center of mass energies cor- 
responding to the 4>, V(3770), or T(AS) resonances, re- 
spectively. The time-dependence of such mesons is com- 
plicated by the issue that the pairs of neutral mesons are 
produced in a coherent wave function consisting of ex- 
actly one |P°) and one |P ) state until one of the mesons 
decays and the correlated wave function collapses. At 
that point in time t\, the second P meson starts to os- 
cillate with mixing frequency AM, until eventually this 
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also decays at some later time t^- The time-difference 
At between these two meson decays replaces the variable 
t used to describe the evolution of uncorrelated mesons. 
The sign of At is taken to be the difference between the 
decay time of a meson into a CP eigenstate minus that 



of the decay into a flavor specific final state (See See- 
Hence events where the CP eigenstate decay is 



tion 



IV) 



the second one to occur have positive values of At, and 
those where the CP eigenstate decay occurs first have 
negative values of At. 



2 T l + |A,|2"- Te U + |A,|* 

h ± Re{\ f ) ' 1 - 12 



The corresponding time-dependence is given by 
r(F ^/) oc e- ri l At l 
T(P°^f) oc e- ri l At l 
where 

h ± = l± e ArAt 
Hence the time-dependent CP asymmetry becomes 

A(At) 



^84fc_ + e—/ 2 ( l^M cos AM At - 2Im ^ 



sin AM At) 



± + R < X 'l h_ + e ArA f /2 cos AMAt + sin AM At ) 



r(Ai)-T(At) _ 2 ArAt/2 (|A/| 2 — 1) cos AM At + 2Im\ f sin AMAi 



r(At) +T(At) 



(1 + |A / | 2 )ft + +2h_Re\ f 



(33) 
(34) 

(35) 

(36) 



and is similar to that for uncorrelated P° production. In 
this case, however, at At = 0, the two P's are completely 
correlated 3 so that the decay of either one is "filtered" 
by the decay mode of the other. When Ar — 0, h + = 2 
and h- = 0. 

For charm decays the measured parameters normally 
used are x and y (or a pair of variables related to x and 
y by a simple rotation), where 



AM , Ar 
, and w = . 

r ' y 2r 



(37) 



Current experimental constraints T give x ~ 0.005 and 
y ~ 0.01. In order to illustrate Eq. (36), the distribution 
.A(At) for D° decays assuming Re\ f = ImX f = l/v2 
is shown in Fig. [2] using x — 0.005 and y = 0.01 [3]. 
It is clear from this illustration that oscillations in the 
charm sector are slow compared with those from Bd or 
B s decays, and the CP asymmetry varies almost linearly 
with At. While an asymmetry is observable, one will re- 
quire large statistics to be accumulated in order to make 
a non-trivial measurement. It should also be noted from 



Eq. ( 36 1 that precise knowledge of both Ar and AM will 



be required in order to translate the slope of the asym- 
metry into a constraint on X f for a given decay channel, 
as indicated by the two curves shown in Fig. [2] 



3 E.g., if the first decays to a CP = — 1 eigenstate then, at At = 0, 
the other has CP = +1 and no odd-CP components will appear 
in its own decay. 



C. Plausibility of measuring time-dependent CP 
asymmetries. 

We will be considering three current or planned 
experimental scenarios where measurements of time- 
dependence of CP asymmetry in charm decays might be 
possible. These correspond to a) LHCb, b) a 2nd genera- 
tion SFF - either SuperB [Ml|| or Belle II [HI [H] run- 
ning e + e~ collisions at the T(4S) or c) Charm threshold 
- Super £? running at the -0(3770) (D°D° threshold where 
the -D's are produced in a coherent state). 

Event yields in each scenario will depend on the specific 
final states considered. Estimates can be made based on 
the proven performance of BABAR and Belle and of the 
current performance of LHCb [33], assuming that the 
current trigger efficiencies can be maintained and that 
the cross section will increase by a factor two for an en- 
ergy increase to 14 TeV from the current 7 TeV. For the 
decays to CP eigenstates (D° -> K + K~ or L>° -> 7T+7T - , 
for example) the yield from a 5 fb _1 sample at LHCb 
is likely to be comparable to that expected from each of 
the SFF's. An upgrade to LHCb could provide a factor 
10 more events. Background levels at LHCb are, how- 
ever, considerably larger than those anticipated at a SFF. 
For modes with higher multiplicity, lower trigger efficien- 
cies will probably contribute to LHCb sample sizes that 
are less competitive with the SFF's, again with larger 
background levels. Event yields at charm threshold will 
be lower, since the SuperB luminosity is expected to be 
smaller at this energy by a factor 10. Backgrounds, how- 
ever, will be lower than at the T{AS). 

To be able to measure t or At, D mesons must be pro- 
duced in flight in the laboratory frame of reference. This 
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FIG. 2: Distribution of A(At) for D° mesons with (top) 
ReXf — ImXf = l/v2 and CP = +1 and (bottom) the ex- 
pected asymmetry for CP = — 1 decay with the same value 
of A/. These distributions assume q/p = 1, and the solid 
(dashed) line corresponds to y — 0.01 (0.00). For Bd decays, 
1.5 full sinusoidal oscillations are observed in the time interval 
presented here. 



means that the flight length of the neutral D mesons un- 
der study needs to exceed the detector resolution associ- 
ated with reconstructing each final state studied. We see, 
however, from Fig. [2] that the asymmetry varies almost 
linearly with decay time so, in principle, asymmetries 
need only be measured in a few regions of At, perhaps 
even for just the two regions At < and for At > 0. 

LHCb should have decay time resolution that is su- 
perior to either of the SFF environments. The D s at 
LHCb are produced with momenta of hundreds of GeV/c 
and the time resolution is generally quite small compared 
to the D° lifetime. The D° mesons at LHCb are pro- 
duced both promptly and in B decays, so care is needed 
to treat each separately. The LHCb trigger has an effi- 
ciency that varies with decay length, and a data driven 
way to measure this variation is required if a system- 
atic limit in precision is to be avoided. Certainly, a finer 
granularity in the time-dependence of any observed CP 
asymmetry is surely possible in the LHCb than in either 
SFF environment. 

Time resolution is more of an issue in the SFF envi- 



ronments. However, prompt D s from e + e~ continuum 
can be cleanly distinguished from those from B decay by 
applying a kinematic cut in momentum. Also, event se- 
lection is not based on decay time, so that efficiency does 
not depend upon decay time. 

D 0, s produced in T(4S) decays have decay lengths 
dominated by the break-up momentum they acquire. As- 
suming the performances of SuperB and Belle II are 
comparable, respectively, to those of BABAR and Belle, 
the decay time resolution for these D°'s are expected 
to be of order one half the D° lifetime. This was suf- 
ficient for both Belle and BABAR to observe mixing in 
the D°-D° system. Therefore, with the yields at the 
Superi? expected to be at least two orders of magnitude 
greater, with backgrounds that are similarly small, and 
with proven data driven techniques to estimate charge 
asymmetry effects, observation of CP asymmetries and 
their time-dependence are at least conceivable. 

At charm threshold, where coherent D° pairs come 
from decays of the "0(3770), the break-up momentum is 
small so that measurements of At rely upon the boost 
from the asymmetric operation of SupcrS. At the T(4S), 
the boost (/?7 ~ 0.23) is such that the decay length (one 
lifetime) for B° mesons is approximately 50/im. Again, 
this corresponds to a time resolution of about one-half 
a B° lifetime. As estimated by the SuperB proponents, 
this is sufficient for measurement of sin 2/3 to a precision 
of 0.1° fI8j. The D° lifetime is, however, 3.8 times shorter 
than that of the B°. At the 0(3770), therefore, a boost 
(Pj) that is approximately four times larger is required 
to maintain the same time resolution - providing the de- 
tector performance is comparable. Decays of D° mesons 
to CP eigenstates have branching fractions about an or- 
der of magnitude larger than those of B° mesons. If the 
larger boost is achievable at the 0(3770), measurements 
of the angle f3 c with a precision similar to that of f3 are 
conceivable. A detailed simulation is required to fully 
understand this, however we discuss results of a simple 



simulation study in Section |VII| focusing on the poten- 
tial for measuring time-dependent asymmetries using the 
D° -> K + K~ and D° -> tt+tt" channels. 



IV. FLAVOR TAGGING 

Flavor tagging is required to synchronize the time t 
in the case of uncorrelated decays and At for correlated 
mesons. Flavor tagging works on the principle of identify- 
ing flavor specific final states that can unambiguously be 
used to determine the flavor of a neutral meson decaying 
into a CP state of interest. A flavor tag has an associ- 
ated probability that the assignment is incorrect. This 
so-called mis-tag probability is denoted by w, and the fig- 
ure of merit used to discuss how useful a particular pro- 
cess or set of channels is for flavor tagging is the dilution 
D = 1 — 2u. More generally one also considers possible 
differences between the mis-tag probability of a particle 
lu and that of the anti-particle uj, where Aw = to — uJ, and 
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the dilution factor becomes D + Auj = 1 — 2uj + Aw. 

An important consideration is that the effect of a non- 
zero value for Auj is an overall shift in CP asymmetries 
at all times and is, therefore, functionally similar to the 
effect of a non-zero value for | X f \ 2 — 1. Uncertainty in Auj 
is, therefore, strongly correlated with that in |A/|. Any 
variation in this quantity with decay time must also be 
well understood if measurements of X f are to be mean- 
ingful. 



A. Flavor tagging of un-correlated mesons 

Flavor tagging of un-correlated D° mesons can be ac- 
complished by identification of "slow" (low momentum) 
pions from the processes D* + — > D°tt + or CP conjugate 
process is D*~ —> D°tt~ . Hence if one can identify a 
sample of events where neutral D mesons originate from 
a D , the charge of the associated pion can be used 
to infer if the D meson is a D° or a D° at the time of 
decay where t = 0. The technical challenge for exper- 
iments is in identifying the so-called bachelor tt^ from 
the D*^ as this has a low momentum and is therefore 
more challenging to reconstruct. At the B factories D* 
from cc continuum can be cleanly separated from D 
from B decay by making a momentum cut above the 
kinematic threshold imposed by the B mass. At LHCb 
the majority of the D* mesons of interest for studying 
CP asymmetries are secondary particles produced in the 
primary decay of a B meson. In either case D* tagged 
events will have non-trivial mis-tag probabilities arising 
from mis-reconstruction, wrongly associated slow pions, 
and from background. A further source of mis-tagging, 
though small, could come from mixing of a D° used for 
tagging. 

One can account for mis-tag probabilities by consider- 
ing the physical decay rates rather than the theoretical 
ones. These are given by 



T Phys^ = (i_ w )r(i)+5jr(i), 

T Phys (t) = w r(t) + (i-aj)r(t), 



(38) 
(39) 



where T(t) and T(t) are from Eqns (|27j and ([28]) . It 
is straightforward to compute the physical CP asymme- 
try by inserting these results into Eq. ( 30 ) . Any preci- 



sion measurement of a CP asymmetry using this method 
would require detailed control of the systematic uncer- 
tainties associated with D* flavor tagging. This provides 
a limit on the ultimate precision attainable for a given 
measurement. 



B. Flavor tagging of correlated mesons 

The set of flavor specific final states of a D meson can 
be used to unambiguously identify if a decay into a CP 



state of interest is that of a D° or a D°. In analogy 
with the methods used for B decay tagging (for exam- 
ple see J7]), one can use a variety of modes for flavor 
tagging D mesons. The advantage of charm over beauty 
can be seen for example in the use of semi-leptonic decays 
for flavor tagging. The decays D — > K^*^i + v account 
for 11% of all D decays, and unambiguously assign the 
the flavor: a D° decay is associated with a £ + in the 
final state, and a D° is associated with a l~ . The cor- 
responding situation for tagging Z?'s from B decays is 
more ambiguous since wrong-sign leptons can arise from 
decays of B's to D* Iv. In addition, the flavour of each 
D° is unambiguously known at At = in the correlated 
case. For uncorrected D 0, s, however, the one decaying 
to a CP eigenstate may have mixed so that its flavour at 
t = is unknown. Thus 11% of all events recorded at the 
■0(3770) can be flavor tagged with a mis-tag probability 
of essentially zero. Events with kaon or pions in the fi- 
nal state can also be used for flavor tagging, however for 
these the mis-tag probability will be non-zero. 

From the perspective of performing a precision mea- 
surement, which will be an inevitable requirement for 
testing the SM, minimization of systematic uncertain- 
ties will be of paramount importance. Here the benefit 
of accumulating data at charm threshold is clear as one 
can choose to restrict the analysis to using only semi- 
leptonic tag decays with an 11% efficiency. In doing so 
an essentially pure CP sample can be reconstructed with 
u) ~ uj ~ 0. 

The viability of including other final states in the tag- 
ging algorithm, for example D° — > K*~ (ir + 1 p + ) etc. in- 
troduces experimental issues that may need to be under- 
stood. These decays can proceed by a tree level Cabibbo 
allowed transition, and the CP conjugate final state can 
proceed via a doubly Cabibbo suppressed transition. 
This introduces an ambiguity in the flavor tag assign- 
ment (hence dilution), and as D mesons can mix there 
are several amplitudes from initial to final state. This 
raises the issue of possible tag-side interference which is 
a well known effect for hadronic B tagging |24) . 

One can account for mis-tag probabilities by consider- 
ing the physical decay rates as a function of At. These 
are given by 



T Phys (At) = (l-E7)r(At)+wF(At), (40) 
ulT(At) + (1 - w)F(At), (41) 



T Phvs (At) 



where T(At) and T(At) are from Eqns (|33j and ((34 1. 
Note that the mistag probabilities are interchanged when 
moving from the uncorrelated (same side tagging) to the 
correlated (opposite side tagging) case. The CP asymme- 
try obtained when allowing for tagging dilution is given 

by 



A _ g^M , (42) 

r v (At)+T ph v s (At) 

(D + Aw)e ArA */ 2 [(|A / | 2 — 1) cos AMAi + 2Im\ f sin AMAt] 

UJ+ / l+ (l + |A / |)V2 + i?e(A / )/ l _ ■ ( } 



Hence a non-zero mistag probability results in a di- 
lution of the amplitude of oscillation, and any particle- 
anti-particle mistag probability differ enc e results in an 
overall offset in the asymmetry 4 . Eq. (43) highlights the 



attraction of using data from charm threshold to min- 
imize systematic uncertainties associated with tagging. 
To a good approximation Alu = 0, and D = 1 for semi- 
leptonic tagged decays, hence the error on A, from this 
source will be relatively small. Furthermore as mentioned 
above, there is only a single amplitude contributing to 
the semi-leptonic tagged side of the event, hence tag-side 
interference is not an issue. Thus if one observes a non- 
zero asymmetry, this can readily be identified as a phys- 
ical effect. For any other tagging category, a significant 
amount of work would need to be done in order to es- 
tablish firstly if the systematic uncertainties were under 
control in terms of tagging performance, and secondly if 
there is a significant issue related to tag-side interference 
that could otherwise manifest large fake signals of CP 
violation. 



V. ANALYSIS OF CP EIGENSTATES 

We have considered a number of two and three body 
CP eigenstate decays of neutral D mesons in order to 
determine the CKM element contributions to the decay 
amplitude, and hence the corresponding weak phase in- 
formation that could be extracted from a given decay. 
The full set of modes is listed in Table [TTJ where we have 
considered contributions from tree, color suppressed tree, 
loop (penguin) and weak-exchange topologies. Possible 
long distance contributions have been neglected in this 
paper. The Feynman diagrams for these topologies, in 
the case of two body final states, are shown in Fig. [3| 

It is clear from Tablellllthat the modes we are consider- 
ing do not contain contributions from all four topologies, 
which simplifies the situation somewhat. One should 
note that the 7r°7r° final state typically consists of four 
photons, however it would be possible to reconstruct a 
vertex and perform a time-dependent analysis for events 
where photon conversion in detector material had oc- 
curred. Also, in about one in 40 instances, one of the 7r°'s 



4 Note that for uncorrelated decays, one interchanges uj and u), 
hence the sign of the Aui terms changes. 



will internally convert in a Dalitz decay 7r° — » e + e~7, in 
which the e-pair with non-zero opening angle will provide 
an excellent location of the vertex position. 

In general we are interested in the value of X f as given 
in Eq. (29) when exploring CP violation. This can be 



written as 



A/ 



(44) 



where 4>mix is the phase of D°D° mixing, and <pcp is the 
overall phase of the D° —> fcp decay, where fcp is a CP 
eigenstate. The amplitude A in general can have contri- 
butions from different topologies, and as a result 4>cp is 
not necessarily directly related to an angle of the charm 
unitarity triangle. This can be seen from the following 



.4 



\T\e v 



\CS\e l 



■ £ \p<\ 

q=d,s,b 



(45) 



where the cf>j , j = T, CS, W, q are phases of the tree, color 
suppressed tree, W exchange and penguin amplitudes re- 
spectively, and the coefficients of the exponentials are the 
magnitudes corresponding to those amplitudes. In gen- 
eral one should note that <pj consists of a strong phase 
(Sj which is invariant under CP) and a weak phase (cj>Y 
which changes sign under CP) , thus cf>j = <pj + Sj . 

If one considers the tree dominated decays such as D — » 
K + K~ , 7r + 7r~, K° K + K~ , and K°tt + tt~ , assuming that 
there is a negligible penguin or color suppressed tree (and 
in the case of a ir + ir~ final state, one also neglects W 
exchange) contribution then it follows that 



A, 



e i<t>Mlx e -2i<p 7 



(46) 



where |T| and the strong phase e %SlT cancel in the ratio of 
A/ A. While this may be adequate for a rudimentary CP 
asymmetry measurement, eventually it would be neces- 
sary to understand the role of the penguin contribution 
to the two body final states, and that of the color sup- 
pressed tree for the three body non-resonant case. It is 
also clear that in order to interpret any CP asymmetry 
measurement in terms of an angle of the charm triangle, 
one needs to obtain a precision measurement of q/p in the 
neutral charm meson system. It should be noted that the 
same arguments also apply for excited states where for 
example pseudoscalar mesons are replaced by vector or 
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TABLE II: CP eigenstate modes considered in this paper indicating the topologies contributing to each process in terms of 
the CKM factors associated with T (tree), CS (color suppressed tree), P q (penguin where q is a down- type quark), and Wex 
(W-exchange) transitions. Blank entries in the table denote that a given topology does not contribute to the total amplitude 
of the decay, and the relative strengths of these amplitudes decrease from left to right. Non-resonant modes are indicated by 
NR in order to differentiate from the resonant contributions with the same final state (but different CP eigenvalue and CKM 
element contribution). 
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axial-vector particles. For final states with two spin one 
particles one must perform an angular analysis in order 
to disentangle CP even and CP odd components of the 
decay. 

In the more general case of two amplitudes contribut- 
ing to the final state (here we consider the case for a tree 
and a single penguin contribution P as a simplification) , 
then 



A, 



re 



e l<PT -\- r ^<PP 



(47) 



where the penguin to tree ratio r = \P\/\T\ is an un- 
known quantity that needs to be evaluated from data. 
If we now return to the amplitudes in Table |TjJ it is 



possible to determine the relative strengths of the differ- 
ent contributions by considering the number of vertices 
in the corresponding Feynman diagram and the CKM 
factors related to these vertices. The distinct products 
of CKM factors appearing in the table are summarized, 
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up to 0(A 6 ), in the following 
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[l-2(p + «?)], (49) 
(50) 



1 
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- ■ - -2 -2 
p — 17] — 7/ — p 



[l-2(p + »?)]. 



(51) 
(52) 



Four of the five amplitudes are complex; V cb V* b has a 
large phase (7J, while V cd V* d 



ub 

and V cd V* s (both have 



the phase of V cd which is (3 C — tt) are related to a small 
weak phase. The remaining term V cs V* d also has a small 
phase, entering at 0(A 6 ) in the amplitude. It is interest- 
ing to note that the amplitude V c bV* b only proceeds via 
a penguin process, and is always accompanied by a tree 
(color allowed or suppressed) and two other penguin am- 
plitudes which will dominate. Hence it is unlikely that 
one will ever be able to collect data with sufficient sta- 
tistical precision to measure j c from processes involving 
a c — > u penguin transition. 

The next most promising phase to measure is asso- 
ciated with transitions mediated by V cd V* d where the 
imaginary component of this amplitude is 0(A 5 ). Modes 
involving this transition at leading order include D — > 
7r + 7r~, p + /0~, h°h°, where h = 7r°,p°,a . These are dis- 
cussed in Sections IV Bl and IV CI 

The combination of CKM elements with the smallest 
phase to this order in A is V cd V* s which is doubly Cabibbo 
suppressed. This CKM factor appears in the W exchange 
amplitudes for D° — > 3K°, however it does so in conjunc- 
tion with other exchange amplitudes that are Cabibbo al- 
lowed. This story is repeated for almost all of the other 
D modes we consider with a neutral or charged kaon in 
the final state. The exceptions D° -> K°f° and K°a° 
have a color suppressed tree proceeding with a CKM fac- 
tor of V cd V* s , and a W exchange amplitude with a factor 
of both V cd V* s and V cs V* d . Hence while the AS ^ 
modes contain weak phase information, it will be dif- 
ficult to experimentally distinguish between the ampli- 
tudes contributing to the decay and extract a precision 
measurement of /3 C . 



A. D° — > K + K and related modes 

D° — > K + K~ measures the phase of V cd V* d only in a 
sub-dominant penguin transition, and is otherwise dom- 
inated by a real tree amplitude with a CKM factor of 
V CS V* S . Hence to first order one would expect to observe 
an asymmetry consistent with the mixing phase 4>mix, 
with no CKM weak phase contribution. This channel 



provides, therefore, a useful cross check of detector re- 
construction and calibration. It also provides measure- 
ments of \q/p\ and 4>mix to complement others that may 
be available. Given that the SM prediction of the asym- 
metry in this channel is small, this is also an ideal mode 
to use when searching for NP. It is interesting to note 
that V cs is complex at C(A 6 ) using the convention of [7]. 
Ultimately a measurement of /3 C could be possible, how- 
ever this is not likely to be the most promising mode to 
measure the angle. 

The same is true for the vector-vector final state 
K* + K*~ . Using the naive factorization framework the 
fraction of longitudinally polarized events Jl in the de- 
cay of a spin zero meson decaying into two vector mesons 
can be estimated as 25 



1 



"i-V 

M 2 ' 



(53) 



where my is the vector meson mass, and M is the mass of 
the decaying parent particle. Using this we can estimate 
f L for D° -> K*+K*- to be ~ 0.77. Hence one would 
be required to perform an angular analysis in order to 
extract CP asymmetry parameters from this decay. 



B. D — > 7r + 7r and related modes 

D° — s- 7r + 7r~ measures the phase of V cd V* d in the lead- 
ing order tree, one of the penguin amplitudes, and the W 
exchange topologies. Of the remaining two penguin am- 
plitudes that contribute to this decay, one is completely 
negligible (mediated by a b quark loop) and the other 
is of the order of A. The non-trivial penguin topologies 
are doubly Cabibbo suppressed loops and proceed at or- 
der A 2 , where as the tree amplitude is singly Cabibbo 
suppressed. A rudimentary measurement of this pro- 
cess could in principle ignore the penguin contribution, 
in which case ImX f ~ sia((f>Mix — 2/3 c )- Thus there will 
be a four- fold ambiguity in any measurement of j3 c . How- 
ever one should note that a more complete analysis would 
be required in order to extract the weak phase and dis- 
entangle the contribution from the c — > s — > u penguin. 

Bigi and Sanda have pointed out [5] that there are 
two Isospin amplitude contributions to D — > Tt + n~ . Ac- 
tually the situation is almost exactly the same as the 
B — > 7T7T, as we have an Isospin 1/2 meson (a B or a D) 
decaying into two pions. The only differences are that, 
in general, we need to assume Ar =^ 0, for charm decays, 
which is a generalization that the existing measurements 
of B° 7T7T have not yet considered, and we neglect the 
W exchange amplitude (which has the same weak phase 
as the tree). The ramification of this is straightforward — 
instead of measuring S and C of Eq. (32 1 in order to de- 



termine the weak phase, one measures the real and imag- 
inary parts of A as given in Eq. (29). One also measures 



7T+A 7T 7r° 



the amplitudes for the Isospin related 7r + 7r 
decays to perform an Isospin amplitude decomposition 
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of 7T7T final states, as described below, in order to disen- 
tangle the phase contribution from the tree and penguin 
amplitudes. 

Similar considerations apply to other final states with 
2-body combinations of n ± , p ± , and 0^(1260). Such 
states that include two spin one particles would require 
an angular analysis in order to disentangle CP even and 
odd parts and correctly measure the time-dependent CP 
asymmetry parameters. For example in the D — > pp 
case, we expect Jl ~ 0.83. As in the B meson system, 
one can apply the same Isospin analysis procedure in or- 
der to bound penguins for D —> pp decays, although one 
should take care to establish whether there is evidence 
of any I — component arising from the finite width of 
the p |26j . Based on the penguin hierarchy observed in B 
decays, we expect that, unless long distance effects play 
an important role in c — > uud transitions, that D — > pp 
might have a smaller penguin contribution than D — > nn. 
If this turns out to be the case, then D — > p + p~ may pro- 
vide a more precise constraint on j5 c than D — > 7r + 7r~, 
and should not be overlooked by experimentalists. It 
should be noted that, while a Quasi- 2-Body approach 
(where the intermediate resonances are treated as par- 
ticles) may be sufficient for a preliminary study, a full 
amplitude analysis would eventually be required in order 
to extract weak phase information from D — > pp decays. 

For decays like, for example, D — > pic, the isospin 
structure can be more complex, in general [57]. We note, 
however, that a complete decay amplitude analysis of the 
7r°7r + 7r~ Dalitz plot has been performed by both CLEO 
|2"5] and by BaBAR [55] and that, in a subsequent isospin 
analysis of this 3-body final state [SO] [31] , it has been 
found that the amplitude is dominated by a single (7 = 0) 
component. This situation is found to be consistent with 
a decay model with no penguin contribution [35] but by 
T, W and CS amplitudes, all with the same phase. This 
makes this channel particularly suitable for extraction of 
f3 c . The BABAR T(4S) sample was very clean and a factor 
five larger than for the tt + it~ channel. A similar state- 
ment can be made for CLEOc running at charm thresh- 
old. For LHCb the trigger is known to be less efficient for 
multi-body final states, which in general produce fewer 
tracks with high transverse momenta to trigger on. As 
a result we do not expect LHCb to be able to make a 
competitive measurement of D — > 7r + 7r~7r° decays when 
compared with the potential of future e + e~ experiments. 
The analysis of this channel is certainly more complex 
than that for 7r + 7r - , but it has been found in both BABAR 
and in Belle experiments that the multi-body channels 
add useful constraints and provide reliable results. 



1. An Isospin analysis of D — > nn and D — > pp decays 

For these decays, the Tree and Penguin decay ampli- 
tudes are distinguished by their isospin changing struc- 
tures. The prescription given here parallels the one de- 
scribed in Ref. 133 which outlines how to measure the 



unitarity triangle angle a from B — ¥ ttt: decays and to 
constrain so-called penguin pollution. Bose symmetry 
dictates that, for either B° or D° decays the two-pion 
final states can be in either an / = or an / = 2 final 
state. In this case, triangular relationships between am- 
plitudes A ij (A^ J ) for D(D) -> tihP decays (h = it or p) 
exist: 



A* 



A+° - A 00 , 



V2 

-±=A- + = A- a -A° a , 
V2 



(54) 
(55) 



where the charges are i,j = +1,-1,0. These two tri- 
angles can be aligned with a common base given by 
A +Q = A , in which case the angle between A^ and 
A + is the shift in the measured phase resulting from 
penguin contributions. 

Obviously, one must measure rates for D° —> h + h~, 
D + -> h + h°, and D° -> h°h° in order to extract the 
weak phase of interest: /3 C . The amplitude of sinu- 
soidal oscillation given in Eq. (301 or (36) is related to 
X f = sin((j)MiX — ZPc,eff)- The proposed Isospin analysis 
would enable one to translate a measurement of Pc,eff to 
a constraint on (3 C , given a precise determination of the 
mixing phase and the amplitudes of D decays to hh final 
states. As final states with more than one neutral par- 
ticle are required for the Isospin analysis, it will only be 
possible to measure the weak phase using D° — > hh de- 
cays in an e + e~ environment. Ultimately the viability of 
this method will depend upon theoretical control of any 
relevant topologies that have been neglected, for instance 
long-distance and isospin-breaking effects. 



D° 



p°p° and related modes 



D° — s- p°p° measures the phase of V c d,V* d via the color 
suppressed tree, one penguin, and W exchange ampli- 
tudes. Of the remaining two penguin amplitudes that 
contribute to this decay, one is completely negligible (me- 
diated by a 6 quark loop) and the other is of the order 
of A. Hence the method to extract the weak phase from 
this decay is a repeat of the situation for D° —> tt + tt~ 
discussed in Section [VB] In order to disentangle the pen- 
guin contribution to the time-dependent CP asymmetry 
measurement, one would have to measure D° — > p + p~, 
which includes two ir° mesons in the final state. So once 
again, this process can only be used to precisely constrain 
the weak phase in an e + e~ environment. It should be 
noted that with p°p°, one can easily measure the time- 
dependent asymmetry, and use the result to reduce the 
number of ambiguities in the D — > pp Isospin analysis. 
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D. New physics 

The topologies summarized in Table [TT] are conve- 
niently categorized in a way where one can envisage dif- 
ferent types of NP affecting the amplitudes contributing 
to the decay rate. NP can manifest itself in any of the 
topologies, and while one normally ignores the possibil- 
ity of NP in tree contributions it is worth noting that 
the measurement of sin 2/3 from B — > J/i/jK are cur- 
rently inconsistent with SM expectations at a level of 
3.2a [5J. This highlights the importance of embarking 
on a quest to measure both the mixing phase and /3 C 
as proposed here. In particular the penguin amplitudes 
could be affected by NP in loop transitions mediated via 
SUSY partners replacing the SM quarks and W ± . Hence 
the modes D° —> h°h a , where h = ir°,p°,(f> are particu- 
larly good candidates to probe NP manifest through this 
mechanism. The remaining modes considered here could 
be used to detect NP contributions from amplitudes that 
compete with the SM tree or exchange amplitudes. In 
general any large observation of CP violation in charm 
decays is expected to be a sign of NP [33]. If one does 
observe a signal, then care must be taken in order to 
disentangle the weak phase of interest from the D° — D° 
mixing phase. This, in turn, will require significantly bet- 
ter measurement of mixing parameters than are currently 
available. 



VI. CONSTRAINING THE SIDES OF THE 
CHARM TRIANGLE 

The charm unitarity triangle given in Eq. (§ can also 
be constrained by measurements of the sides, essentially 
magnitudes of the elements of the CKM matrix. The dif- 
ference in the lengths of the two long sides V* d V c d and 
V* S V CS must be able to accommodate the geometry of 
the third side indicated in Fig. [l] While the direct mea- 
surement of CP violating effects is the focus of this pa- 
per, the indirect measurements required to constrain the 
shape of the triangle independently of CP asymmetry 
measurements are also important and worthy of a men- 
tion. We briefly examine each of these elements in turn 
in the following to highlight how one can increase current 
knowledge of the triangle via indirect measurements. 

| V u d | : This has been precisely measured using nuclear 
beta decay and the experimental level of precision 
reached is at the level of 0.022% [55] . 

| V us | : This quantity can be measured precisely in kaon 
decays, however that has reached a natural conclu- 
sion of being dominated by systematic uncertain- 
ties. The level of precision reached for this quan- 
tity by averaging results from kaon and r decays 
is 1% [35]. Future precision measurements of \V US \ 
may be possible via studies of r decays into final 
states with charged kaons. Thus the SuperB and 



Belle II experiments will be required to improve our 
knowledge of this quantity. 

|V^b|: The limiting factor for improving constraints on 
this element comes from a combination of theoret- 
ical and experimental issues relating to B decays 
into semi-leptonic final states related to b — > u tran- 
sitions. While there has been a lot of work in this 
area, there is still a lot of room for improvement 
both in terms of theoretical and experimental de- 
velopments. The current level of uncertainty ob- 
tained for this quantity is 11% [35]. From the ex- 
perimental perspective the inclusive and exclusive 
results obtained for \ V u t,\ are not in good agreement 
with each other [3B] . Thus the SuperB and Belle II 
experiments will be required to improve our knowl- 
edge of this quantity. 

iV'cdh Precision measurements of semi-leptonic D de- 
cays can improve our knowledge of | V c d\ beyond the 
current level of precision (4.8% [35]). This measure- 
ment can be improved upon by the BES III exper- 
iment at IHEP, and also by the Superi? and Belle 
II experiments. SuperB will have the advantage of 
being able to accumulate at data sample fifty times 
larger than BES III at charm threshold. It is un- 
likely that Belle II would ultimately be competitive 
with a measurement of \ V c d\ as that experiment has 
no plans to run at charm threshold. 

|V^ S |: The most precise determinations of | | come 
from measurements of semi-leptonic D s decays. 
The current level of precision obtained for |V CS | is 
3.5%. This can be improved by the BES III experi- 
ment at IHEP, and also by the SuperB experiment , 
using data collected just above charm threshold. 

|V c b|: The limiting factor for improving constraints on 
this element comes from a combination of theoret- 
ical and experimental issues relating to B decays 
into semi-leptonic final states related to b — > c tran- 
sitions. While there has been a lot of work in this 
area, there is still a lot of room for improvement 
both in terms of theoretical and experimental de- 
velopments. The current level of precision achieved 
by measurements of \V c b\ is 3.2% 35 . From the ex- 
perimental perspective the inclusive and exclusive 
results obtained for \ V c b\ are not in good agreement 
with each other [36] . Thus Super £? and Belle II will 
be required to improve our knowledge of this quan- 
tity. 

\V u d\ is the most precisely constrained quantity re- 
quired to reconstruct the triangle using the sides hav- 
ing been measured to 0.022%. Hence improved measure- 
ments of this quantity will not play an important role 
in improving our understanding of the charm triangle. 
All of the other quantities are known to precisions of the 
order of 1 — 10%. Thus in order to improve indirect con- 
straints of the charm triangle, (i) we need to wait for the 
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SuperB and Belle II experiments to improve the limiting 
factors in terms of measuring the above quantities, and 
(ii) the corresponding theoretical developments should 
also be pursued in order for experiment to remain a lim- 
iting factor. It should also be noted that the BES III 
experiment will be able to improve the precision of mea- 
surements of | V c d | and | V cs | from semi-leptonic D and D s 
decays before the SFF's start collecting data. 

Interestingly enough the quantities |I4f>| and \ V c b\ also 
currently limit the precision of the sides constraint of 
the unitarity triangle for B decays, and again the only 
routes to experimental improvements on that test are via 
SuperB and Belle II. 



VII. NUMERICAL ANALYSIS 

In this section we compare the three experimental sce- 
narios, (i) charm threshold (ii) the SFF's at T(45), and 
(hi) LHCb, relating to the measurement of CP violation 
in D° —> fcp decays, where fcp is a CP eigenstate. We 
neglect resolution effects related to the reconstruction of 
vertices in the detector and translation of this spatial 
distance into values of At or t. Finally, based on the ex- 
pectations from these simulations, we discuss the direct 



constraint on the apex of the cu triangle in Section VII E 



For the numerical analysis and the extrapolation to 
the expected precision in /3 c , e / / , we generate a set of one 
hundred Monte Carlo data samples in each experimen- 
tal scenario, For Super_B running at charm threshold we 
do this for both semi-leptonic and also kaon decays as 
tags. In each sample, we generate D° and D° events 
with no time-integrated asymmetry, each according to 
their respective time dependences described in Sec. III. 
We simulate effects of mis-tagging either D° or D , then 
perform a binned fit to the resulting asymmetry given in 
Eq. (43) and the corresponding form in terms of t. In 



these fits, arg(A/) and |A/| are allowed to vary and the 
values for w and Aw are fixed at those used in the event 
generation. We repeat this analysis for different possible 
values of the phase arg(A/) from —10° to +10° in 10° 
steps. As a figure of merit for each experimental sce- 
nario, we take the average uncertainty, cr^, in this phase 
from the 100 fits, observing that this is consistent with 
the spread of central values from the individual fits. 



A. Charm threshold 

D meson pairs produced at the 0(3770) are quantum- 
correlated, so that the time evolution is given by 
Eqns (p33"|) and (34). If one accounts for tagging dilu- 



tion, then the time-dependent CP asymmetry is given 
by Eq. (43). On restricting time-dependent analyses to 
the use only of semi-leptonic tagged decays, the asym- 
metry simplifies as there is no dilution, since both w and 
Aw terms can be neglected, and any systematic uncer- 
tainty in the asymmetry arising from D ~ 1 becomes 



small. Furthermore the e + e~ — ¥ -0(3770) environment 
is extremely clean, so that systematic uncertainties from 
background contributions is also small and under control. 
These are important points to stress as we know that the 
CP phase of interest is expected to be small, hence in 
order to make a precision measurement the systematic 
uncertainties must be minimized. 

With 500 fb -1 of data at charm threshold one can 
expect to accumulate approximately 1.8xl0 9 D meson 
pairs. With a data sample of 281pb~ 1 CLEO-c obtain 
89 D° —> 7r + 7r~ candidates with the other D meson de- 
caying semi-leptonically into X + ev e . Their efficiency 
for such events is 50% |37J . Assuming the same effi- 
ciency applies 5 , we anticipate that SuperB could record 
158, 000 Xev e tagged D° — > tt + tt~ events, correspond- 
ing to 489500 events when using the full set of K(*kve 
tagged events, I = e, /i. We expect about three times 
the number of events for D° K + K~ . Figure [4] shows 
the results obtained for the average uncertainties in the 
phase arg(A^) = <f> = 4>mix + 4>cp as a function of that 
phase. 
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FIG. 4: The uncertainty in the measured phase <j> = 4>mix — 
2j3 c as a function of the value of <j> for (squares) D° — > -k + tt~ 
(triangles) D° — > K + K~ decays at charm threshold with 
500 ft^ 1 of data, assuming that the mis-tag probability is 
negligible, and only using the full set of semi-leptonic tagged 
decays. 



These results are only for semi-leptonic tags. We also 
consider use of hadronically tagged events, for example 
D° —> K~X (K + X), where X is anything, which corre- 
spond to 54% (3%) of all neutral D meson decays. From 
these modes alone, one would expect w ~ 0.03, and that 
the asymmetry in particle identification of K + and K~ in 
the detector will naturally lead to a small, but non-zero 
value of Aw. We expect that there would be approxi- 
mately 2.2 million kaon tagged D° — > tt + tt~ events in 
5001b- 1 at charm threshold. Using these data alone, one 



Preliminary studies indicate that this is a reasonable assumption. 
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would be able to measure ^ to a precision of 4°. Hence, 
if one combines the results from semi-leptonic and kaon 
tagged events, a precision of er^ ~ 3.4° is achievable. 
This represents a significant improvement in precision 
over just using semi-leptonic tagged events 6 . 



B. Uncorrelated decays at the T(45) 

The scenario at the T(4S) is somewhat more com- 
plicated than the situation encountered at the ^(3770). 
Firstly, in order to remove background from D mesons 
produced in B meson decay, one restricts the analysis to 
mesons with high momentum. In addition to non-trivial 
backgrounds, one also has to consider non-zero tagging 
dilution, where the asymmetry is similar to that given 
in Eq. (43), but with t substituted for At, and u) and 57 



interchanged (hence a sign flip for the Aw terms). Thus 
it is not obvious that Aw can be neglected, and indeed 
D/l. BABAR recorded 30,679 D* tagged D° -> Tr+Tr" 
events at the T(4S) in 384 fb -1 of data [35], with a pu- 
rity of 98%, and where the mis-tag probability for these 
events is ~ 1% [39 . From this we estimate that one could 
reconstruct 6.6xl0 6 D* tagged D° — > h + it~ events in a 
data sample of 75 ab" 1 . We obtain the sensitivities for 
arg(A / ) = <f> as a function of the phase shown in Fig. H 
assuming this yield and dilution. 



C. Uncorrelated decays at LHCb 

The final scenario considered is that of measuring time- 
dependent asymmetries from uncorrelated D mesons in 
a hadronic environment. Preliminary time-integrated re- 
sults from CDF [40] and LHCb [23] indicate that such 
a measurement is possible. Dilution and background ef- 
fects will, however, be more severe in this hadronic envi- 
ronment than at an e + e~ machine. The measurement of 
| X f | is expected to be dominated by such systematic un- 
certainties, though arg(A f ) may be less affected, provided 
that any variation of u> or Auj as a function of decay time 
can be carefully controlled. It is not clear at this point 
what the ultimate precision obtained from LHCb will be. 
The best way to ascertain this would be to perform the 
measurement. 

Based on the result in Ref. [23] we estimate that LHCb 
will collect 7.8xl0 6 D* tagged D° -> 7r+7r~ decays in 
5 fb -1 of data, based on an initial 37pb _1 of data. Based 
on the data shown in the reference, we estimate a purity 
of ~ 90% and ui ~ 6%. From these values, we obtain the 
sensitivities for arg(A/) = (/> as a function of the phase 
shown in Fig. [6] assuming this yield and mis-tag proba- 
bility. 
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FIG. 6: The uncertainty in the measured phase cj> = 4>mix ~ 
2/3 c as a function of the value of <f> for D° —¥ tt + tv~ decays at 
LHCb with 5fb _1 of data, assuming cj = uj — 0.06. 



FIG. 5: The uncertainty in the measured phase cj> — 4>mix — 
2/3 c as a function of the value of <\> for D° — > tt + ty~ decays at 
the T(4S) with 75 ab -1 of data assuming to = u> — 0.01. 

D. Summary of sensitivity estimates 



To offset the aforementioned background and dilution 
issues, the increased boost in this T(4S) scenario does 
slightly reduce the effects of time resolution that are ig- 
nored in our analysis here. 



6 Use of the K tag events will introduce tag side interference. For 
Bd analyses this amounts to a few parts per mille, but it will 
need to be evaluated for the specific D modes that are used. 



Measurements of arg(A / ) in all scenarios will require 
good knowledge of the D°D° mixing parameters. These 
should be available from SFF's running at T(4S) and 
from LHCb. Superi?, for example, expects [18] to mea- 
sure these with precisions of a few times 10~ 4 for x and y, 
~ 1.5° for 4>mix and a few % in \q/p\. More information 
comes from the D° — > K + K~ sample, and it is likely that 
LHCb will further improve these parameters. Hence, in 
a fit combining these measurements, it will be possible 
to separate out contributions from the mixing and weak 
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phase in D a —> ir + Tr~ decays. More accurately, the dif- 
ference between phases of X f measured in D° — > K + K~ 
and D° — > ir + ir~ decays is <ficp = ~^Pceff- If 1°°P con- 
tributions can be well measured and both long-distance 
and weak exchange contributions are negligible, then this 
constraint can be translated into a measurement of (3 C . 

In order to relate the measurement of the weak phase, 
(3 c ,eff of X f to p c , one needs to measure a set of Isospin 
related D — > hh decays. This is something that will re- 
quire the e + e~ environment, as it will not be possible for 
LHCb to reconstruct D° -> 7r°7r°, or D° -> p + p~. Both 
D + — » 7r + 7r° and D + — > p + p° would also be challenging 
measurements for LHCb. Nonetheless a search for CP 
violation in D° — » ir + ir~ and D° — > p°p° decays and a 
measurement of 4>mix at LHCb would be of considerable 
interest. 



The corresponding sensitivity estimates for arg(A / ) in 
the different scenarios considered are summarized in Ta- 
ble [Hi] It should be recalled, however, that effects from 
time resolution, or from time-dependencies in efficiency 
or mis-tag rates are neglected here. We estimate that it 
should be possible to measure 4>cp to ~ 2.6° using this 
approach. Assuming that penguin contributions can be 
measured precisely, then the error on /3 C that could be 
obtained by Super B would be ~ 1.3°. LHCb will require 
input from SuperB on the decay modes with neutral par- 
ticles in the final state in order to translate a measure- 
ment of Pc,ejf to one on j3 c . Further work is required 
to understand how penguins and other suppressed am- 
plitudes affect the translation of P c ,eff to /3 C , however it 
is clear that there will be a significant contribution from 
penguins given the size of the D° —> 7r°7r° branching frac- 
tion. It is worth noting that BABAR and Belle could be 
able to make a measurement of /? c ,e// with a precision 
of ~ 25°, using the nominal values of x and y measured 
for charm mixing available today. We have highlighted 
several decays that could be used to measure this angle, 
including D — ► tttt, pn, pp, aiir, K°f°, and K°a°. Ulti- 
mately it will be important to measure /? c ,e// in each of 
these modes in order to cross-check the consistency of all 
of the measurements, constrain NP, and bound possible 
corrections to the CKM mechanism. 



TABLE III: Summary of expected uncertainties from 500 fb 1 
of data at charm threshold, 75 ab _1 of data at the T(iS), and 
5 fir 1 of data from LHCb for D° -> K+ K~ and D° -> ty+-k~ 
decays. The column marked SL corresponds to semi-leptonic 
tagged events, and the column SL+K corresponds to semi- 
leptonic and kaon tagged events at charm threshold. 



Parameter 


SL 


SuperB 
SL + K 


T(4S) 


LHCb 


4>(-k-k) = arg(A,r T ) 


8.0° 


3.4° 


2.2° 


2.3° 


4>(KK) = &rg{X KK ) 


4.8° 


2.1° 


1.3° 


1.4° 


4>CP = 4>KK — (pTTTT 


9.4° 


3.9° 


2.6° 


2.7° 


Pc,eff 


4.7° 


2.0° 


1.3° 


1.4° 



Here we have concentrated on the determination of the 
phase arg(Aj), and one should not neglect the fact that 
we are also able to constrain \X f \ using these same mea- 
surements. An observation of |A/| ^ 1 in data would con- 
stitute the measurement of direct CP violation in a given 
decay channel. We estimate that it should be possible to 
measure |A/| with a statistical uncertainty of 1 — 4% at 
the future experiments discussed above, though we note 
that this is limited by any uncertainty in Aw. This is 
smallest in SuperB running at charm threshold, but is 
statistically limited less in other scenarios. 

If one compares the relative power of data from charm 
threshold with that from the T(AS), it is clear from Ta- 
ble^ 
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that 75 (50) ab _1 of data at the T(4S) is equivalent 
to approximately 1.2 (0.8) ab" 1 at charm threshold. It 
is interesting to note that SupcrB proponents expect to 
accumulate 500 fb -1 of data at charm threshold in only 
three months, whereas 75 ab _1 would require five years of 
running at nominal luminosity. The time-scale involved 
for the Belle II experimental run at the T(AS) is similar 
to the SuperB one. 



E. Constraint on the cu triangle 

It is possible to constrain the apex of the cu triangle 
in Figure [T] by constraining two internal angles, or by 
measuring the sides. If one considers the representation 
where the baseline V* S V CS is normalized to unity, then 
the angles at vertices corresponding to the coordinates 
(0,0) and (1,0) are f3 c and a c , respectively. The con- 
straint on the apex of the cu triangle can be obtained 
using the CKM prediction of 7 C = (68.4±0.1)° (from the 
Bd triangle), and any future measurement of j3 c . The 
7 C constraint is essentially a straight line in the complex 
plane containing the cu triangle. As is the case with the 
Bd triangle, there are multiple solutions for the apex of 
the triangle. Even a rudimentary constraint on /3 C , made 
by establishing that /3 c ,eff is compatible with zero, would 
constitute a test of the SM. A precision measurement of 
Pc,eff would require a detailed treatment of theoretical 
uncertainties to determine if any small deviation from the 
expected value of /3 C was due to new physics, or compat- 
ible with the SM. This is an area that will require work 
in the future. We are currently working on determining 
the effect of penguin pollution in D — > hh decays. In 
addition to this effect, other potential sources of theo- 
retical uncertainty that may be relevant include isospin 
breaking effects, long distance topologies or failure of the 
factorization hypothesis. The coordinates of the apex of 
the triangle are given by 



X + iY 



A 2 X 5 (p + if]) 



X - A 3 /2 — A 5 (l/8 + A 2 /2) '' 



(56) 



neglecting contributions from all higher orders in A. 
Given that the apex of the bd triangle is p + if], one 
can over constrain the SM by testing the prediction of 
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X + iY from existing constraints on the apex of the bd 
triangles. We find that 

X = 1.00025, (57) 
Y = 0.00062, (58) 

using the existing constraints on the Wolfenstein param- 
eters. 

In order to measure (3 c ,eff one needs to precisely con- 
strain <Pmix- The current method to measure the mix- 
ing phase is via a time-dependent Dalitz Plot analysis 
of D decays to self conjugate final states. Here we pro- 
pose to use a time-dependent analysis of decays such as 
D — » K + K~, which have an overall phase dominated by 
the mixing phase in the SM assuming the CKM param- 
eterization, and rate larger than the tttt channel. Hav- 
ing determined 4>mix one can then decouple the mixing 
phase contribution in D — > ir + n^ decays, and by per- 
forming an Isospin analysis one can translate a measure- 
ment of X f into a constraint on f3 c ,eff- Alternatively one 
can use a model independent measurement of the mixing 
phase, to decouple (f>Mix and /3 c , e // from the measure- 
ment of X f . One would have to control both theoretical 
and systematic uncertainties to below one per mille in 
order to be sure of measuring a non-zero value of /3 C . At 
this time it is unclear if this will be achievable, however 
the Superf? experiment has the added advantage of be- 
ing able to study the time-dependence in two ways, and 
hence may be able to avoid limitations inherent to the 
D* tagged analyses. 



VIII. B d DECAYS 

The effect of a non-zero Ar on the time-dependent CP 
asymmetry distribution is an alteration of the phase of 
oscillation, and of the amplitude of the oscillation as a 
function of t or At. Until now all time-dependent CP 
asymmetry measurements in B^ decays have assumed 
that Ar = 0, which was a reasonable assumption based 
on theoretical expectations. However it should be noted 
that it is possible to bound the systematic uncertainty in 
the measurement of the unitarity triangle angles a and /3 
by making this assumption using the known experimental 
constraint on sign(ReX f )AT/T = 0.010 ± 0.037 [35]. If 
one compares the asymmetry obtained assuming Ar cor- 
responding to the experimental bound, then it is possible 
to estimate the bias and systematic uncertainty time- 
dependent asymmetry measurements made in B decays 
arising from the assumption that Ar = 0. 

We have performed a Monte Carlo based simulation 
for the scenario of S = 0.7 and C = 0.0 taking the uncer- 
tainty in Ar to be Gaussian. The ratio of amplitudes for 
the first maximum/minimum obtained as an estimate of 
the systematic effect on S = sin 2/3 is 0.007 ± 0.027, and 
the corresponding distribution is shown in Pig. [7] This is 
comparable to the statistical uncertainty in sin 2/3 mea- 
surements [ini [T7] . 



Moving onto the measurements related to a, if one 
considers B° — > tt + it~ decays, where S — —0.65 ± 0.07 
and C = -0.38±0.06 gHHS], then the systematic uncer- 
tainty in the measurement of S and C is 0.009±0.032. In 
this case, the systematic effect resulting from the assump- 
tion that Ar = is also non-trivial, but does not domi- 
nate the total uncertainty. The most important channel 
for the constraint on a is however B° — > p + p~ where S = 
-0.05 ± 0.17 and C = -0.06 ± 0.13 [H US]. The corre- 
sponding systematic effect on S and C is —0.008 ±0.038, 
which is currently small compared to the experimental 
determination of those quantities. Therefore, while the 
Ar = bound may impact upon the /3 constraint im- 
posed on the unitarity triangle, it will have little effect 
on the measurement of a. 

Therefore current and future experiments aimed at 
performing a precision measurement of time-dependent 
CP asymmetries should also strive to increase the preci- 
sion of the bound on Ar to ensure that this systematic 
effect docs not dominate future measurements. 




Amplitude Ratio 

FIG. 7: The bias on S = sin 2/3 obtained from a Monte Carlo 
based simulation resulting from the assumption that Ar = 0. 
The amplitude ratio plotted is that of the maximum time- 
dependent amplitude accounting for a non-zero AT to that 
where Ar = 0. 



IX. B a DECAYS 

Oscillations in B s decays are extremely fast relative to 
Bd and D mesons, and so neither SuperB or Belle II are 
expected to be able to perform time-dependent asymme- 
try measurements in B s decays. It should however be 
noted that if these experiments were to accumulate large 
samples of events at the T(5S), then the distribution of 
events as a function of At would contain information on 
both the real and imaginary parts of X f . Hence some 
information from CP asymmetries related to the time- 
dependent measurements being done at hadron collider 
experiments would be measurable in an e + e~ environ- 
ment. This was also discussed in [33] in the context of 
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measurements of B s — > J/ip<fi at SuperB. This issue 
is particularly relevant for final states including neutral 
particles such as B s — > rf<j>, the B s equivalent to the 
most precisely measured golden b — > s penguin mode 
B° — > rj'K°. It would be extremely challenging to study 
this mode a hadronic environment and so the best way 
to study CP violation in this mode would be using data 
collected at the T(5S). 

Other interesting decays to study are B s — > 
pK°,DfK T , and B s -!• D<\> as these measure 7 [44H46] . 
It would be interesting to compare the values obtained 
from a B s decay with the result from the Bd — > DK ap- 
proach currently being used by experiments. It should 
be noted that LHCb should be able to perform time- 
dependent measurements of these modes. Finally, as 
noted in Ref. [35], the channel B s — > tt q K2. is equiva- 
lent to the channel Bd — > 7r + 7r~. Therefore it would 
be interesting to attempt to measure X f for this decay. 
Given the ir° in the final state, and lack of information 
to constrain a primary vertex, this could be an excellent 
candidate for SuperB or Belle II to study. 

X. CONCLUSIONS 

We have outlined the formalism required to exper- 
imentally measure time-dependent CP asymmetries in 
charm decays using correlated D°D° decays as well as 
D° mesons tagged from D* decays, and discussed the 
benefits of studying a number of different CP eigenstates. 
The important points to note are that one can use K + K~ 
decays to measure the mixing phase quite precisely and 
other decays can be used to constrain the angle (3 c ,eff 
which is related to the cu unitarity triangle. These ob- 
servables are also sensitive to possible enhancements from 
new physics. A data sample of 500 ftU 1 collected at 
charm threshold would provide a sufficient test to con- 
strain any potential large NP effects. Similar measure- 
ments would also be possible using D* tagged decays at 
SuperB, Belle II and LHCb. From event yields currently 
available, we expect the statistical precision in the mea- 
sured phase at SuperB to be slightly better than results 
from a 5 fb _1 LHCb run. As the cu and the bd unitar- 
ity triangles are related, the measurements proposed here 
provide a new set of consistency checks on the unitarity of 
the CKM matrix that can be performed using D decays. 
Measurements of the sides of these triangles would enable 
a further, indirect cross-check on the validity of this ma- 



trix. Only the SuperB experiment will be able to make 
a complete set of the measurements required to perform 
direct and indirect constraints of both triangles. As f3 c is 
an extremely small angle, its determination will be lim- 
ited by theoretical and systematic uncertainties. Superi? 
has a potential advantage over other experiments as it 
will be able to collect data at charm threshold with a 
boosted center of mass, as well as being able to explore 
effects using neutral mesons from D* tagged events. Data 
from charm threshold will be almost pure, with a mis-tag 
probability of ~ for scmi-lcptonic tagged events, which 
could be advantageous if systematic uncertainties dom- 
inate measurements from T(45 l ) data and from LHCb. 
The ultimate theoretical uncertainty in relating /? c ,e/ / to 
j3 c needs to be evaluated. A measurement of |A / | 7^ 1 
could also signify direct CPV . 

We also point out that precision measurements of time- 
dependent asymmetries in Bd decays require improve- 
ments in our knowledge of ATs d - The current exper- 
imental constraint on this observable translates into a 
systematic effect of the order of 0.007 ± 0.027, which is 
comparable with the current experimental sensitivity on 
sin 2/3 from BABAR and Belle. We have also computed 
the systematic effect of assuming Ar s d = for measure- 
ments of a from B° — > ir + ir~ and B° — > p + p~ decays, 
which is negligible for existing measurements. 

It may be possible to measure the real and imaginary 
parts of X f from a simplified time-dependent analysis of 
B s decays at SuperB and Belle II without the need to 
observe oscillations. While the approach outlined would 
not be competitive with modes that could be measured 
in a hadronic environment, it would provide unique ac- 
cess to observable channels that would be inaccessible to 
the Tevatron and LHCb. The prime example is that of 
B s — > rj'tfi, which is the direct analog of the most pre- 
cisely measured B® — > s penguin mode B® rj' K° from 
the B factories. 
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